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Abstract The extended nebulae formed as pulsar winds expand into their surround- 
ings provide information about the composition of the winds, the injection history 
from the host pulsar, and the material into which the nebulae are expanding. Obser- 
vations from across the electromagnetic spectrum provide constraints on the evo- 
lution of the nebulae, the density and composition of the surrounding ejecta, the 
geometry of the central engines, and the long-term fate of the energetic particles 
produced in these systems. Such observations reveal the presence of jets and wind 
termination shocks, time-varying compact emission structures, shocked supernova 
ejecta, and newly formed dust. Here I provide a broad overview of the structure of 
pulsar wind nebulae, with specific examples from observations extending from the 
radio band to very-high-energy y-rays that demonstrate our ability to constrain the 
history and ultimate fate of the energy released in the spin-down of young pulsars. 



1 Introduction 

The basic structure of a pulsar wind nebula is determined by the spin-down energy 
injected by the central pulsar and the interaction of the nebula with the interior re- 
gions of the supernova remnant (SNR) in which it evolves. Losses from synchrotron 
radiation in the nebular magnetic field, whose strength depends both on the nature 
of the injected wind and on the evolving size of the PWN, inverse-Compton (IC) 
scattering of ambient photons by the energetic electron population within the neb- 
ula, and adiabatic expansion as the nebula sweeps up the surrounding supernova 
ejecta, all combine to determine the emission structure and long-term evolution of 
the nebula. (See (18 1 for a review.) Multiwavelength observations of PWNe provide 
crucial information on the underlying particle spectrum and strongly constrain both 
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the magnetic field strength and the stage of evolution. Of particular interest is the 
spectrum of low-energy particles contained in the PWN. These retain the history of 
early energy losses as well as possible signatures of features in the pulsar injection 
spectrum. 

Complex structure in PWN spectra can originate in a number of ways that are as- 
sociated with the long-term evolution as well l32l[T9l . In particular, when the reverse 
shock compresses the PWN, the resulting increase in the magnetic field initiates a 
new epoch of fast synchrotron cooling that adds to the population of low energy 
particles. Models which do not account for the reverse shock interaction will thus 
underestimate this population and its associated radiation. In addition, recent studies 
of the spectrum immediately downstream of the wind termination shock show that 
the injection spectrum itself can deviate significantly from a simple power law ||38ll , 
and particle-in-cell simulations of the acceleration process produce a Maxwellian 
population with a power law tail l40l . Any such structure in the injected particle 
spectrum imprints itself on the broadband emission of the entire nebula. The re- 
sulting breaks or regions of curvature in the PWN spectrum, and the frequencies at 
which they are observed, depend upon the energy at which features appear in the 
electron spectrum as well as the means by which the photons are produced (e.g. 
synchrotron radiation or IC emission). To fully understand the nature of the particle 
injection, as well as the long-term evolution of PWNe, it is thus crucial to study the 
emission structure over the entire electromagnetic spectrum. 



2 Dynamical Evolution of PWNe 

The evolution of a PWN within the confines of its host SNR is determined by both 
the rate at which energy is injected by the pulsar and by the density structure of 
the ejecta material into which the nebula expands. The location of the pulsar itself, 
relative to the SNR center, depends upon any motion given to the pulsar in the form 
of a kick velocity during the explosion, as well as on the density distribution of the 
ambient medium into which the SNR expands. At the earliest times, the SNR blast 
wave expands freely at a speed of ~ (5 — 10) x 10 3 km s~ 1 , much higher than typical 
pulsar velocities of ~ 400 — 500 km s . As a result, for young systems the pulsar 
will always be located near the SNR center. 

The energetic pulsar wind is injected into the SNR interior, forming a high- 
pressure bubble that expands supersonically into the surrounding ejecta, forming 
a shock. The input luminosity is generally assumed to have the form (e.g. Il30l ) 




jn+l) 
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Fig. 1 Composite image of Crab Nebula. X-ray emission from Chandra is concentrated in a 
central jet/torus structure (shown in blue), while optical emission from HST (shown in yellow and 
red), and the infrared emission from Spitzer (shown in purple) dominates the exterior regions where 
the nebula sweeps up ejecta material. 



is the initial spin-down time scale of the pulsar. Here Eq is the initial spin-down 
power, Pq and Po are the initial spin period and its time derivative, and n is the so- 
called "braking index" of the pulsar (n = 3 for magnetic dipole spin-down). The 
pulsar has roughly constant energy output until a time To, beyond which the output 
declines fairly rapidly with time. In the spherically symmetric case, the radius of the 
PWN evolves as 

r^^lse^e^m;;^ 5 ( 3) 

where E$n is the energy released in the explosion and M e j is the mass of the ejecta 
iflOl . Thus, at least at early times when the pulsar input is high, the PWN expan- 
sion velocity increases with time. The sound speed in the relativistic fluid within 
the PWN is sufficiently high (c s = c/y/3) that any pressure variations experienced 
during the expansion are quickly balanced within the bubble; at early stages, we 
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thus expect the pulsar to be located at the center of the PWN. The pressure bal- 
ance within the PWN results in a termination shock where the energetic pulsar wind 
meets the more slowly-expanding PWN. Studies of this region of the PWN provide 
the most direct information available on particles that are being freshly injected into 
the nebula. 

The geometry of the pulsar system results in an axisymmetric wind |25:|, form- 
ing a torus-like structure in the equatorial plane, along with collimated jets along 
the rotation axis. The higher magnetization at low latitudes confines the expansion 
here to a higher degree, resulting in an elongated shape along the pulsar spin axis 
for the large-scale nebula Bl l44ll . This structure is evident in Figure Q] where X- 
ray and optical observations of the Crab Nebula clearly reveal the jet/torus structure 
surrounded by the elongated wind nebula bounded by filaments of swept-up ejecta. 
A thorough review of the MHD-driven jet/torus structure is presented by N. Buc- 
ciantini in these Proceedings. 

Beyond the PWN, the outer blast wave of the SNR drives a shock into the sur- 
rounding interstellar/circumstellar medium (ISM/CSM), forming a shell of hot gas 
and compressed magnetic field. As the shell sweeps up additional mass, and decel- 
erates, a reverse shock (RS) propagates back into the expanding ejecta. When the 
SNR has swept up an amount of mass that is roughly equal to the mass of the ejecta, 
the evolution approaches the Sedov phase in which the radius is described by 



where «o is the density of the ambient medium. This is illustrated in Figure [2] where 
dashed lines show Rpwn and Rsnr as a function of time using Equations 3 and 4. 
Here we have assumed M e j = 5M Q and E$n = 10 51 erg, and have used a range 
of values for Eq and uq. As indicated, eventually the PWN would overtake the SNR 
boundary under such conditions. However, because the initial expansion of the SNR 
occurs more rapidly than in the Sedov phase, and the injection rate from the pulsar is 
highest at the earliest times, the actual behavior differs from that shown. In addition, 
as the reverse shock of the SNR propagates inward, it eventually reaches the PWN 
boundary and begins crushing the nebula. This is indicated in the solid curves in 
Figure [31 where the early -phase behavior as well as the reverse shock propagation 
is taken into account (see |fl9l ). Ideally the PWN radius decreases until the nebula 
pressure is sufficiently high for the system to rebound. 

The crushing of the PWN results in an increase in the magnetic field. This causes 
a rapid burn-off of the most energetic particles in the nebula. The PWN/RS inter- 
face is Rayleigh-Taylor (R-T) unstable, and is subject to the formation of filamentary 
structure where the dense ejecta material is mixed into the relativistic fluid |6|. The 
nebula subsequently re-forms as the pulsar injects fresh particles into its surround- 
ings, but a significant relic nebula of mixed ejecta and relativistic gas will persist. 

In cases where the pulsar and its PWN have moved considerably from the center 
of the PWN, or in which the SNR has evolved in a medium of nonuniform den- 
sity, the reverse shock will interact with the PWN asymmetrically, encountering one 
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Fig. 2 Time evolution of the SNR and PWN radii for a range of values for the ambient density and 
initial spin-down power of the pulsar. The solid curves correspond to models from Gelfand et al. 
(2009) fl9l using En = 10 40 erg s~\ M ej ■ = 8M @ , n Q = 0.1 cm- 3 ,and£ 51 = 1. See text description 
for details. (From |39|. Reproduced by permission of the AAS.) 



portion of the nebula well before another. This results in a complex interaction that 
leaves a highly distorted relic nebula that may be highly displaced from the pulsar 
position (6|. 

In the later stages of evolution, the PWN can expand to a very large size, with 
a correspondingly lower magnetic field strength. As described below, such PWNe 
may be best identified through their high-energy y-ray emission, with only weak 
X-ray emission observed close to the pulsars. 



3 Spectral Evolution of PWNe 

As particles are injected from a pulsar into its PWN, the resulting emission is de- 
termined by the evolved particle spectrum and magnetic field, as well as the energy 
density of the ambient photon fieldQ The injected spectrum is often characterized 
as a power law: 

Q(E e ,t)=Q (t)(E e /E )- a (5) 
The integrated energy in the electron spectrum is then 



1 In addition, emission from ejecta gas and dust swept up by the expanding nebula can be significant 
- and even dominant in some wavebands. 
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Fig. 3 Synchrotron (left) and IC (right) emission (for scattering off of the CMB) from a PWN at 
ages of 1000 (solid), 2000 (dotted), and 5000 (dashed) years. Here we have assumed £51 = 1, M e j = 
8M ( t), and «o = 0.1 cm~ 3 for the SNR evolution, and n = 3, £0 = 10 40 erg s~', and To = 500 yr for 
the pulsar. For the wind, we assume that 99.9% of the energy is in the form of electrons/positrons 
with a power law spectrum with 7 = 1.6. 



where O" is the ratio of the spin-down power injected in the form of Poynting flux to 
that in the form of particles. It is important to note one expects the postshock flow 
of particles to be characterized by a Maxwellian distribution accompanied by a non- 
thermal tail - a result confirmed by recent particle-in-cell simulations of relativistic 
shocks ll40l . The development of the nonthermal tail depends on the shock condi- 
tions, and the total residence time of the particles in the acceleration region, which 
may vary at different locations in the inner nebula. Thus, it is reasonable to expect 
that the injected spectrum may actually deviate considerably from a pure power law 
form, a point that recent observations may be beginning to illustrate, as we discuss 
below. 

The resulting emission spectrum is found by integrating the electron spectrum 
over the emissivity function for synchrotron and IC radiation using, respectively, 
the nebular magnetic field and spectral density of the ambient photon field. As il- 
lustrated in Figure [3] the build-up of particles in the nebula results in an IC spec- 
trum that increases with time. The synchrotron flux decreases with time due to the 
steadily decreasing magnetic field strength associated with the adiabatic expansion 
of the PWN. This behavior is reversed upon arrival of the SNR reverse shock (not 
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shown in Figure), following which the nebula is compressed and the magnetic field 
strength increases dramatically, inducing an episode of rapid synchrotron losses. 
Upon re-expanding, however, IC emission again begins to increase relative to the 
synchrotron emission. At the latest phases of evolution, when the nebula is very 
large and the magnetic field is low, the IC emission can provide the most easily- 
detected signature. As described below, such behavior is seen for a number of PWNe 
that have been identified based on their emission at TeV energies, and for which only 
faint synchrotron emission near the associated pulsars is seen in the X-ray band. 

The broadband spectrum of a PWN, along with the associated dynamical infor- 
mation provided by measurements of the pulsar spin properties, and the size of the 
PWN and its SNR, place very strong constraints on its evolution and on the spectrum 
of the particles injected from the pulsar. Combined with estimates of the swept-up 
ejecta mass, this information can be used to probe the properties of the progenitor 
star and to predict the long-term fate of the energetic particles in the nebula. Re- 
cent multiwavelength studies of PWNe, combined with modeling efforts of their 
evolution and spectra, have provided unique insights into several of these areas. 



4 Case Studies 

A complete summary of recent results on PWNe is well beyond the scope of this 
paper. Below we discuss four distinct systems for which recent multiwavelength 
observations have begun to probe both the detailed structure and the underlying 
particle population in these nebulae. 



4.1 3C58 

3C 58 is a flat-spectrum radio nebula (a ps 0.1, where S v « V~ a ) for which upper 
limits based on IRAS observations indicate a spectral break between the radio and 
infrared bands J2J]. The PWN has often been associated with SN 1 181 EH; the low 
break frequency would then suggest an extremely large magnetic field (> 2.5mG) 
if interpreted as the result of synchrotron break. This has resulted in a number of 
different interpretations, including the suggestion that the pulsar in 3C 58 under- 
went a rapid decline in its output at some early epoch |2T1 , the possibility that the 
low-frequency break is inherent in the injection spectrum from the pulsar, and the 
suggestion that 3C 58 is not actually associated with SN 1181, but is an older nebula 

una. 

The X-ray emission from 3C 58 (Figure|4|i is dominated by a power law compo- 
nent, typical of synchrotron emission. However, a faint thermal component is clearly 
detected in the outer regions of the PWN, and also contributes to the interior regions 
|8]|37]|22]. The temperature is ~ 0.25 keV, and enhanced abundances of Ne and Mg 
are observed, indicative of ~ 0.5Mq of ejecta that has been swept up by the PWN. 
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Fig. 4 Top: Chandra image of 3C 58. The pulsar is at the center, and is surrounded by a elongated 
compact nebula with a curved jet extending to the west. A softening of the spectrum with radius is 
observed - an effect resulting from both synchrotron aging of the electrons and the presence of a 
soft thermal shell. Bottom: Chandra image of the pulsar in 3C58, and its associated torus and jet 
(left) and the IRAC 8 (jm image of the same region. The elliptical region indicates the torus, and 
has the same center and size in each image. 



This is much larger than the expected mass if the PWN is associated with SN 1181, 
and suggests a larger age for the system [ 1 1 1. 

To investigate the evolution of 3C 58, we have carried out Spitzer observations 
using IRAC 11381 . These observations reveal the PWN in both the 3.6 and 4.5 fim 
bands, representing the first detection of synchrotron emission from this important 
young PWN anywhere in the five decades of frequency separating the radio and 
X-ray bands. The morphology of the IR emission from 3C 58 is strikingly similar 
to that seen in the radio and X-ray bands. The emission extends all the way to the 
radio boundaries, indicating that no synchrotron loss breaks occur below the band, 
and some regions of enhanced or diminished emission match well with those seen 
in the other bands (notably the large cavity on the eastern side), suggesting that we 
are observing primarily synchrotron radiation. Optical filaments in 3C 58 l37l[33ll . 
which presumably originate from supernova ejecta overtaken by the expansion of 
the PWN, do not show a good spatial correspondence with the radio or IR struc- 
tures, suggesting that the IR emission is not dominated by dust contributions. This 
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Fig. 5 The flux of all of 3C 58 (upper) and its torus (lower), plotted from the radio to the X-ray 
band. While the torus is not detected in the radio band, the IR data require a flattening of the X- 
ray spectrum when extrapolated back to the longer wavelength band. (From |38| . Reproduced by 
permission of the AAS.) 

is similar to results from Spitzer observations of the Crab Nebula B2| . where emis- 
sion in the IRAC band is also identified primarily with synchrotron radiation. 

The IRAC data also reveal emission from the torus surrounding the pulsar in 
3C 58 |36) in all four bands (FigureB)0 Optical emission from the torus is detected 
as well l35ll . These observations provide new constraints on the evolution of the 
particles as they flow from the termination shock in 3C 58. There is little question 
that this emission is synchrotron in nature; there is insufficient dust in the environ- 
ment of the pulsar termination shock to provide a shocked dust component to the 
emission. The spectrum (Figure|5]l requires a break of some sort between the IR and 



2 We also show preliminary results from MIPS observations at 24 fim where we also detect the 
torus. 
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Fig. 6 ROSAT PSPC image of Vela SNR. The single contour represents the outer boundary of the 
radio nebula Vela X. The point source at the northern extreme of Vela X is the Vela pulsar and the 
box inside the PWN indicates the XMM MOS region of the cocoon shown in the inset. 



X-ray bands, suggesting that the synchrotron loss break appears just above the IR 
band. Most importantly, these results indicate that the spectrum of particles injected 
into the PWN through the termination shock does not follow an unbroken power 
law. As a result, structure in the PWN spectrum is, at least in part, the imprint of 
structure from the injection region. 

From Figure [5] it is clear that additional observations of 3C 58, and particularly 
its torus, at longer wavelengths will be crucially important to understand the nature 
of the injected particles and the subsequent long-term evolution of the PWN. Deep 
observations of the central regions of other PWNe are clearly of importance as well. 



4.2 VelaX 



Located at a distance of only 280 pc, the Vela SNR houses a young pulsar that pow- 
ers the extended nebula Vela X. This nebula lies within a limb-brightened shell of 
thermal X-rays. The outer shell is cool, and the exceptionally low foreground ab- 
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Fig. 7 XMM-Newton spectrum from the "cocoon" region in Vela X. The best-fit model, shown in 
black, is composed of two components - a thermal plasma with enhanced, ejecta-like abundances 
(light dashed curve) and a power law (dark dashed curve). (From [23]. Reproduced by permission 
of the AAS.) 

sorption (A^ ~ l-5 x 10 20 cm~ 2 ) allows us to see strong emission lines from O, 
Ne, and Mg. In X-rays, the overall brightness asymmetry of Vela is evident (Fig- 
ure[6]l. The SNR is much brighter in the northeastern hemisphere, toward the Galac- 
tic plane. This is apparently the result of large-scale inhomogeneities in the ISM, 
with «o ~ 0.06 cm~ 3 in the south Q, and «o ~ 1 cm~ 3 on the north side of the SNR 

m. 

Radio observations of the PWN fl28l reveal a morphology concentrated to the 
south of the pulsar itself, suggesting that the nebula has been disrupted by the im- 
pact of the reverse shock which propagated more rapidly from the northeast due to 
the higher ambient density in this direction. Higher resolution radio images show a 
network of filamentary structure in the PWN H151 . possibly formed by R-T instabil- 
ities in this interaction with the reverse shock. ROSAT observations of the Vela X 
region l27l reveal a large emission structure - the so-called "cocoon" - extending 
to the south of the pulsar. The region is characterized by a hard spectrum and ap- 
pears to lie along a bright elongated radio structure. ASCA observations established 
a two-component X-ray spectrum with the hard component adequately described by 
either a power law or a hot thermal plasma (29]. The PWN is observed at energies 
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Fig. 8 Broadband spectral model consisting of synchrotron emission in the radio and X-ray bands 
accompanied by IC emission in the VHE y-ray band. The upper panel shows a model with two 
spectral breaks in the electron spectrum. Models with a single break (lower panel) either underpre- 
dict the radio emission, or produce a radio flux whose spectral index does not agree with observa- 
tions. (From [231 . Reproduced by permission of the AAS.) 



up to ~ 200 keV with BeppoSAX l26l . and observations with H.E.S.S. |2) reveal ex- 
tended VHE y-ray emission with the brightest emission concentrated directly along 
the cocoon. 

Our initial studies of a small region along the cocoon |[23l reveal a bright X-ray 
structure shown as an inset to Figure [6] The emission is concentrated into several 
distinct regions, at least some of which appear to be filamentary structures. The 
integrated emission from these regions is characterized by two distinct components 
- a power law with a spectral index of ~ 2.2 and a thermal plasma with enhanced 
abundances of O, Ne, and Mg, presumably associated with ejecta that has been 
mixed into the PWN upon its interaction with the reverse shock (Figure |7J. 

It is of particular interest that our broadband modeling of the nonthermal emis- 
sion from this central region of Vela X indicates a disconnect between the radio- 
emitting particles and those that produce the X-ray and TeV y-ray emission (23). A 
model with a single spectral break between the bands either underpredicts the radio 
flux or produces an incorrect radio spectral index (Figure [8]), while a model with 
two breaks can satisfy the data. Treating the poorly-characterized low-energy elec- 
tron component as a separate population of particles, de Jager, Slane, & LaMassa 
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Ifl2ll showed that enhanced IC emission from this component could be expected 
in the GeV band. This has now been confirmed with observations by AGILE l3T1 
and Fermi JT|. Interestingly, a similar excess is observed in Fermi observations of 
HESS J1640-465 (39J, another evolved PWN that appears to have undergone a RS 
interaction (see Section 4.4). 

The nature of the low-energy particle spectrum in Vela X is poorly understood, 
but the Fermi studies in particular suggest a difference in the cocoon emission from 
that of its surroundings; the LAT emission appears to be concentrated distinctly to 
the west of the TeV emission. Whether or not this emission component is somehow 
associated with the reverse shock interaction is not clear. An XMM Large Project 
to map a significant portion of Vela X is underway to investigate the distribution of 
nonthermal particles and thermal ejecta in the nebula. 



4.3 G327.1-1.1 

G327. 1 — 1 . 1 is a composite SNR with a bright central PWN whose structure is com- 
plex in both the radio and X-ray bands. As shown in Figure|9l the SNR is character- 
ized by a faint radio shell accompanied by a bright radio PWN. The PWN is offset 
from the SNR center, and a finger-like structure extends toward the west/northwest. 
The morphology is suggestive of a PWN that has been disrupted by an asymmet- 
ric reverse shock interaction that has arrived preferentially from the northwestern 
direction. 

Detailed Chandra observations 11431 appear to confirm this scenario in detail! A 
compact X-ray source resides at the tip of the radio finger, and a trail of nonthermal 
X-ray emission extends from the source back into the radio nebula. The compact 
X-ray source itself is resolved, with possible evidence of a jet or torus structure, and 
the source is embedded in a cometary nebula whose structure is suggestive of a bow 
shock. A pair of prong-like structures originate from the vicinity of the compact core 
and extend out to the west/northwest. Their axes are not aligned with the compact 
core, and they do not appear to be jets from the pulsar. The most unusual feature 
in G327. 1-1.1 is a large bubble-like structure that extends out from the prongs. The 
bubble is very faint, but it is apparent in images taken with both Chandra and XMM. 

The spectral index in the tail of X-ray emission that extends into the brighter radio 
PWN appears to vary from ~ 1.8 to ~ 2.1, although significantly better statistics 
are needed to map this fully. Faint loop-like structures, possibly associated with 
Rayleigh-Taylor filaments or magnetic loops, are evident in the southern part of the 
tail and they appear to extend from radio structures inside the relic PWN. Deeper 
X-ray observations are required to better characterize these and, more importantly, 
to study the detailed properties of the complex structures seen in the immediate 
vicinity of the pulsar, and to investigate the apparent upstream diffusion of particles 
that form the bubble-like structure. 



3 See also the contribution by T. Temim in these Proceedings. 
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Fig. 9 Composite radio (red and orange) and X-ray (blue) image of G327. 1 — 1.1. The outer radio 
shell defines the SNR boundary, while the bright central nebula is the PWN. The compact source 
to the northwest of the nebula is the neutron star. It is embedded in a cometary-shaped structure 
accompanied by a tail of X-ray emission extending into the radio nebula, as well as prong-like 
structures (seen in inset) that appear to inflate a faint bubble in the northwest. 



The overall X-ray morphology of the PWN in G327. 1-1.1 presents several chal- 
lenges. It appears clear from the displacement of the radio nebula that the PWN has 
undergone an interaction with the SNR reverse shock, and that this shock arrived 
earlier from the northwest, possibly as a result of pulsar motion in this direction. 
In such a scenario, we expect the ongoing pulsar wind production to begin forming 
a new PWN around the pulsar, shaped by the surrounding pressure conditions. Ra- 
dio emission from the displaced relic nebula will persist, but the compression from 
the reverse shock will temporarily increase the magnetic field in the nebula, causing 
rapid synchrotron losses for the more energetic particles. The X-ray emission is thus 
expected to be concentrated closer to the pulsar, as observed. Spectral steeping of the 
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X-ray spectrum in the direction of the relic nebula is expected if the synchrotron loss 
timescale is shorter than the particle flow timescale; mapping this spectral evolution 
along the extended X-ray tail thus constrains the conditions in the relic nebula. In 
the direction from which the reverse shock propagated, the structure of the medium 
is complicated. The density in the immediate post-shock region is enhanced, but 
declines downstream due to adiabatic expansion of the SNR. At least qualitatively, 
this could provide the environment in which freshly-injected wind from the pulsar 
inflates the observed bubble. 

The cometary feature surrounding the putative pulsar complicates the above pic- 
ture. The morphology is suggestive of a bow shock that forms when the pulsar mo- 
tion exceeds the sound speed in the the ambient medium, yielding a structure quite 
different from that of a static PWN. The wind termination shock in such systems 
is compressed in the direction of motion, and extended in the backward direction. 
X-ray emission is then observed between the termination shock and the contact dis- 
continuity, and forms three distinct structures: the "head" that surrounds the pulsar; 
an enhanced region behind the pulsar, associated with the termination shock; and 
an elongated tail where the swept-back wind is concentrated lfT7l l9l. For G327.1- 
1.1, such a geometry would require a significant pulsar velocity component per- 
pendicular to the plane of the sky, but the overall velocity that would be required 
(770 km s ) is not unreasonable. 

Given the similarity between G327. 1—1.1 and Vela, in the context of both being 
systems in which an interaction between the PWN and the SNR reverse shock has 
occurred, it is of considerable interest to determine whether or not G327. 1 — 1.1, 
like Vela X, produces y-rays. Preliminary investigation of Fermi-LAT data indicates 
faint emission that is positionally coincident with G327. 1 — 1.1, but further analysis 
is required to assess this in detail. 



4.4 HESS J 1640 -465 

HESS 11640—465 (see Figure [TUb is an extended source of very-high-energy y-ray 
emission discovered with the High Energy Stereoscopic System (H.E.S.S.) during a 
survey of the Galactic plane 0. Centered within the radio SNR G338. 3-0.0 l34l 
the deconvolved TeV image of the source has an RMS width of 2.7 ±0.5 arcmin 
ifTBTl . HI measurements show absorption against G338.3— 0.0 out to velocities cor- 
responding to the tangent point, indicating a distance of at least 8 kpc 1241 . and 
thus implying a rather large size for the PWN (RpwN > 6.4<tfio pc, where c/jo is the 
distance in units of 10 kpc). X-ray observations with XMM |fl6l and Chandra l24l 
establish the presence of an accompanying X-ray nebula and an X-ray point source 
that appears to be the associated neutron star. The point source is offset from the 
center of the PWN, and both are offset from the center of the SNR, suggesting that 
an asymmetric interaction with the SNR reverse shock has occurred. 

We have investigated Fermi-LAT data acquired from the region surrounding 
HESS 11640—465 and detect the source with high significance [39]. The spectrum 
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Fig. 10 Fermi-LKT image of HESS J1640— 465. The solid (cyan) circle indicates the uncertainty 
in the centroid of the Fermi-LAI source, the magenta dashed circle indicates the 95% encircled flux 
distribution of the H.E.S.S. image, and the white dashed circle indicates that for 3EG J1639— 4702. 
The white contours outline radio emission from G338.3— 0.0 while the black contours at the cen- 
ter outline extended X-ray emission observed with XMM. A compact X-ray source detected with 
Chandra resides within the X-ray contours. (From |39| . Reproduced by permission of the AAS.) 

is well-described by a power law with r = 2.30 ±0.09 and a F(> 100 MeV) = 
2.8 x 10~ 7 photons cm~ 2 s _1 . We have modeled the emission assuming a 1-zone 
model in which particles are injected into the nebula with a simple power law dis- 
tribution. We use a radius of Rsnr ~ 11.6dio P c f° r G338.3— 0.0, based on radio 
observations. The observed extent of HESS J1640— 465 constrains the radius of the 
PWN to Rpwn > 6.4ifio pc. As indicated in Figure|2] where a horizontal line indi- 
cates the radius of G338.3— 0.0, reasonable values for the SNR and PWN parameters 
indicate that the SNR reverse shock has almost certainly begun to impact the PWN. 

The broadband emission model results are shown in Figure QT| where we plot 
the Fermi and H.E.S.S spectra along with the radio upper limit from GMRT obser- 
vations 11201 and spectral confidence bands derived from Chandra 1241 . The black 
curves represent the model prediction for the synchrotron (left) and IC (right) emis- 
sion that best describes the X-ray and Te V y-ray spectra, similar to results from l24l ; 
the parameters for the model are summarized in the caption. As seen in Figure [TT] 
this model significantly underpredicts the observed Fermi-LAI emission. Our spec- 




Fig. 11 Electron spectrum (upper) and broadband emission model (lower) for HESS J1640— 465 
assuming the evolutionary history described in the text. The black curves represent a PWN with 
an age T = 10 kyr, and B(T) = 5/iG, assuming £ — 4 x 10 36 erg s and an injection spectrum 
with a = 10~ 3 , 7= 2.5, and E m ; n = 115 GeV. The light curves represent the scenario with a low- 
energy Maxwellian electron component replacing the low-energy portion of the electron power- 
law spectrum. The mean temperature for the IR and optical photon fields are 15 K and 5000 K, 
respectively, and the energy densities relative to the CMB are 4 and 1.15. The dashed curve in the 
upper panel represents the truncated portion of the power law that was replaced by a Maxwellian. 
The dashed curve in the lower panel represents a model for which all of the y-ray emission results 
from pion decay. (From | 39|. Reproduced by permission of the AAS.) 



tral fits can formally accommodate up to about ~ 20% of the observed flux in a 
pulsar-like component characterized by a power law with an exponential cutoff en- 
ergy between 1 and 8 GeV, and there are several known radio pulsars located within 
the error circle of HESS J1640— 465 that could potentially produce observable y-ray 
emission. Even in this case, however, the Fermi-LKY emission still greatly exceeds 
the predicted flux from HESS J 1640-465. 
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As described in Section 4.2, simple power-law models for the particles in Vela 
X, another evolved PWN, fail to reproduce the observed broadband spectrum. The 
presence of an excess population of low-energy electrons is inferred, and models 
for the IC scattering of photons from this population predict an excess of y-rays in 
the GeV range. Motivated by these results from Vela X, we modified the evolved 
power law spectrum from our model for HESS J1640— 465 by truncating the lower 
end of the power law and adding a distinct low-energy component. Based on results 
from simulations of shock acceleration ll40l . we chose a Maxwellian distribution 
for this population. Our resulting (ad hoc) particle spectrum is shown in the upper 
panel in Figure QT| and the resulting broadband emission is shown in the lighter 
(magenta) curves in the lower panel. Here we have adjusted the normalization of 
the Maxwellian to reproduce the emission in the Fermi-LAT band, which is pro- 
duced primarily by upscattered infrared (IR) photons from local dust. We find a 
mean value of y « 2 x 10 5 for the electrons in the Maxwellian component, and 
roughly 5% of the total electron energy in the power law tail, consistent with results 
from particle-in-cell simulations. Recent work by Fang & Zhang (2010) |fl4| uses 
a similar input distribution to successfully model the emission for several PWNe 
including HESS J1640-465. However, their results for HESS J1640-465 under- 
predict the observed GeV emission from this source, apparently due to use of a 
slightly lower bulk Lorentz factor and a larger fraction of the total energy in the 
power law tail than we have used in this analysis. 

An alternative scenario for the y-ray emission is that it arises from the SNR 
itself, and not the PWN. The dashed blue curve in Figure QT| represents a model 
for the emission from the collision of protons accelerated in the SNR with ambient 
material, leading to y-rays from the production and subsequent decay of neutral 
pions. Assuming a shock compression ratio of 4 and that 25% of the total supernova 
energy appears in the form of relativistic protons, an ambient density no ~ 100 cm -3 
is required to produce the model shown in the Figure. This is much higher than can 
be accommodated for the observed size of the SNR and the lack of observed thermal 
X-ray emission from the SNR. Such high densities are found in dense molecular 
clouds, suggesting that the y-rays could be produced by particles that stream away 
to interact with high-density material outside the SNR. However, only the most 
energetic particles can escape the acceleration region, which is in conflict with the 
proton spectrum we require to match the data. Moreover, the observed Te V emission 
appears to originate from within the SNR boundaries, making such an escaping- 
particle scenario appear problematic. Based on this, along with the lack of a spectral 
cutoff that might suggest emission from a central pulsar, we conclude that the GeV 
y-ray emission most likely arises from the PWN. 



5 Summary 



The broadband spectra of PWNe provide information about both the structure and 
evolution of these objects. New multiwavelength observations have begun to probe 
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PWNe from the sites of particle injection to the ejecta-laden outer boundaries, pro- 
viding crucial input for modeling these systems. Observations in the y-ray band have 
uncovered previously unknown systems in the late phase of evolution, while X-ray 
observations continue to provide detailed information about the geometry and the 
composition of the pulsar winds. These observations continue to inform theoretical 
models of relativistic shocks which, in turn, have broad importance across the realm 
of high-energy astrophysics. At the same time, these recent results have pointed the 
way to new and deeper observations of PWNe across the electromagnetic spectrum. 
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